Odd-particle systems in the shell model Monte Carlo: circumventing a sign problem 
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We introduce a novel method within the shell model Monte Carlo approach to calculate the 
ground-state energy of a finite-size system with an odd number of particles by using the asymptotic 
behavior of the imaginary-time single-particle Green's functions. The method circumvents the sign 
problem that originates from the projection on an odd number of particles and has hampered direct 
application of the shell model Monte Carlo method to odd-particle systems. We apply this method 
to calculate pairing gaps of nuclei in the iron region. Our results are in good agreement with 
experimental pairing gaps. 
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Introduction. The shell model Monte Carlo (SMMC) 
approach [lHj] has been used successfully to calculate sta- 
tistical properties of nuclei [5|-|7| within the framework of 
the configuration-interaction shell model. Recently, this 
method has also been applied to trapped cold atom sys- 
tems @, Q ■ The SMMC method enables calculations in 
model spaces that are many orders of magnitude larger 
than those that can be treated by conventional diagonal- 
ization methods. 

For typical effective nuclear interactions, the SMMC 
method breaks down at low temperatures because of the 
so-called fermionic sign problem, leading to large statis- 
tical errors. In the grand-canonical ensemble the sign 
problem can be avoided by constructing good-sign inter- 
actions that include the dominant collective components 
of effective nuclear interactions (loj . The remaining part 
of the effective interaction can be accounted for by using 
the method of Ref. 0. 

In finite-size systems, such as nuclei, it is necessary 
to use the canonical ensemble, in which the number of 
particles is fixed. This particle-number projection gives 
rise to an additional sign problem when the number of 
particles is odd, leading to a rapid growth of statistical 
errors at low temperatures even for good-sign interac- 
tions. Consequently, it has been a major challenge to 
make accurate estimates for the ground-state energy of 
odd-particle systems in SMMC. Accurate ground-state 
energies are necessary for the calculation of level densi- 
ties and pairing gaps (i.e., odd-even staggering of binding 
energies). 

Here we develop a method based on the asymptotic 
behavior of the imaginary-time single-particle Green's 
functions of an even-particle system to calculate ground- 
state energies of neighboring odd-particle systems. This 
method is somewhat similar in spirit to a technique used 
in lattice quantum chromodynamics to extract hadron 
masses (see, e.g., in Ref. [Tlh . We apply our Green's func- 
tion method to calculate pairing gaps of nuclei in the iron 
region using the complete fp + 59/2 shell model space. 

Green's functions in SMMC. The SMMC method is 
based on the Hubbard-Stratonovich representation of the 



imaginary-time propagator, e~P H = J D[a]G(o-)U a (/3), 
where /J is the inverse temperature, H is the Hamilto- 
nian, D[a\ is the integration measure, G(a) is a Gaussian 
weight, and U a ((3) is the propagator of non- interacting 
nucleons moving in external auxiliary fields a that de- 
pend on the imaginary time r (0 < r < /3). The canonical 
thermal expectation value of an observable O is given by 
(6) = J D[a]G(a)Tr A [dU a (0)]/ J Dio-lG^TrAU^P), 
where Tr_4 denotes a trace over the subspace of a fixed 
number of particles A. In actual calculations we project 
on both proton number Z and neutron number N, and 
in the following A will denote (Z, N). 

For a quantity X a that depends on the auxiliary fields 
cr, we define 

= p[4WA m 

a ~ jD[a]\W(a)\^ ' ( > 

where W{a) = G{a)Tv A U a and $ ff = W{a)/\W[a)\ is 
the sign. With this definition, the above thermal expec- 
tation of an observable O can be written as (O) = (O) , 
where (6) a = Tr A [6u a (f3)}/Tr A U a (/3). In SMMC we 
choose M samples according to the weight function 
|W(it)|, and estimate the average quantity in (JT|) by 

For an even number of particles with a good-sign in- 
teraction, the average value of the sign $ CT remains close 
to 1. However, when the number of particles is odd, the 
average sign decays towards zero as the temperature is 
lowered. This leads to rapidly growing errors, hampering 
the direct application of SMMC at low temperatures for 
odd-particle systems. 

For a rotationally invariant and time-independent 
Hamiltonian, we define the following scalar imaginary- 
time Green's functions 

GAt) = Tr^ e-m ' ( 2 ) 

where v = (nlj) labels the nucleon single-particle orbital 
with radial quantum number n, orbital angular momen- 
tum I and total spin j. Here T denotes time ordering 
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and a um (r) 



e TH a v 



-tH 



is an annihilation operator 



of a nucleon at imaginary time t (— ft < t < f3) in & 
single-particle state with orbital v and magnetic quan- 
tum number m (— j < m < j). 

Using the Hubbard-Stratonovich transformation, the 
Green's functions defined in ([5]) can be written in a form 
suitable for SMMC calculations 



G v {t) 



£[U CT (T)(I-(p) CT U !( , m forr>0 
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for r < 
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where we have used the notation in Eq. ([T]). Here U CT (r) 
and I are matrices in the single-particle space repre- 
senting the propagator U a {r) and the identity, respec- 
tively. {p) a is a matrix in the single-particle space whose 
vm^v'm! matrix element {p V m.v'm')a is defined in terms 



of the one-body density operator p 
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Assuming A is an even-even nucleus, A± = (Z, N ± I ) 
are neighboring odd-even nuclei with odd number of neu- 
trons. We denote by Jn (n = 0, f , 2, . . .) the n-th excited 
state with total spin J and define the energy differences 
AEj(A±) = E J0 (A±) - E 00 (A), where E,, n {A) is the 
energy of the state Jn with particle number A. Assum- 
ing that the ground state of the even-even nucleus has 
spin zero, AEj(A±) is the energy difference between the 
lowest state of a given spin J in the odd-even nucleus 
A± and the ground state of the even-even nucleus A. 
Assuming that the ground state of the odd-even nucleus 
A± is J — j, where j is one of the single-particle or- 
bital spin values, its corresponding energy is given by 
Eg S (A±) = E Q0 {A) + AE min (A±), where AE min is the 
minimum of AEj(A±) over the possible values of j. 



The neutron Green's function G v {t) that corresponds 
to an orbital with angular momentum j can be written 
as 



G v (t) = C(P)e 



-AEj(A±)\r\ 
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where the + (— ) subscript should be used for r > 
(r < 0) and C(/3) is a r- independent con- 
stant. Rj™ (A± , v) are scaled weights defined by 
Ril l '(A + ,v) = |(jV||at||Jn)| 2 /|(j0||4||00)| 2 and 
#fn '(■*-; i/) - |(JV|| a ,||Jn)| 2 /|(jO||a,||00)| 2 , where 
(J'n'||aJ,|| Jn) and {J'n'Wa^Jn) are reduced matrix ele- 
ments of a\, and a v between the state Jn in A and the 
states J'n' in A+ and A-, respectively. 

When all terms in the summation on the r.h.s. of Eq. 
(0| are small, the Green's function can be well approxi- 
mated by a single exponential, G v {t) ~ e ~ AEj ( A± ^' r \ . In 
this asymptotic regime for r, we can calculate AEj(A±), 
and hence E gs (A± ) from the slope of In G v (r) . This is 
the method we use here to calculate the ground-state en- 
ergy of odd-A nuclei with odd number of neutrons. The 
ground-state energy of odd-A nuclei with odd number 
of protons can be similarly calculated using the proton 
Green's functions. 

In principle, the asymptotic regime is accessed in the 
limit (3 — > oo. However, in a shell-model Hamiltonian 
with discrete, well separated energy levels, only a few 
transitions give significant contributions. If the relative 
contribution from the sum in Eq. Q is less than a few 
percent, then (assuming that r| ~ I MeV) the sensitiv- 
ity of the slope of In G„ (r) to this contribution is about 



a few tens of keV, which is comparable to our target ac- 
curacy. For low- and medium-mass nuclei, we expect the 
energy differences to be > 1 MeV and the scaled weights 
to be much smaller than one. Thus, calculations with (3 
of a few MeV -1 and with an asymptotic regime of r ~ 1 
MeV should be sufficient. This can be validated explic- 
itly in sci-shell nuclei (see below) , whose Hamiltonian can 
be diagonalized numerically. For larger model spaces, it 
is not possible to calculate explicitly the corrections in 
the sum of Eq. Q , and the asymptotic region has to be 
determined by the goodness of the linear fits to In G v (r). 

Results. We first tested the Green's function method 
in s<i-shell nuclei and then applied it to medium-mass 
nuclei in the complete (pf + gg/2) shell. In these nuclei, 
we carried out calculations for several values of j3 in the 
range 3 MeV -1 < (3 < 4 MeV -1 . For each /3, we calcu- 
lated G v {t) for a range of values of r in steps of 1/32 
MeV -1 . We chose the asymptotic region in t such the 
linear fits to In G v {t) have a \ 2 P er degree of freedom 
~ I or less in all cases considered. We find that a good 
asymptotic region is 0.5 MeV -1 < r < 2 MeV -1 . 

Within the asymptotic region, we fit a straight line to 
lnG„ for each possible subset of points in t for which 
G v (t) has been calculated. The mean and standard de- 
viation of the slopes so obtained are used to estimate 
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FIG. 1: The absolute value of logarithm of the Green's func- 
tion ([2]) for the neutron orbital v — lpz/2 in 56 Fe (lower curve, 
t > 0) and 58 Fe (upper curve, r < 0) at f) = 4 MeV -1 . The 
solid blue lines are linear fits for 0.5 MeV -1 < \t\ < 2 MeV -1 . 
The insets show the deviations from these linear fits. 



AE m i n (A±) and its statistical error, respectively, at each 
j3. A weighted average of the results at different values 
of /3 is then taken. 

In a few selected cases, we also performed calculations 
for larger values of (3 (i.e., (3 > 4 MeV -1 ), and found the 
corresponding values of AE m i n (A±) to be consistent with 
those obtained in the region 3 MeV -1 < p < 4 MeV -1 . 
This indicates that for the model spaces and particle 
numbers considered, the above chosen values of (3 are 
sufficiently large to isolate the ground state of the corre- 
sponding even-even nucleus. 

For a given odd system (an odd-even nucleus) there are 
two neighboring even systems (even-even nuclei) , and our 
method can be used by starting from either of the even 
systems. Unless noted otherwise, the results we report 
here are the average of both of these calculations. 

To test the validity and accuracy of our method, we 
performed calculations in the sd shell using a schematic 
good-sign Hamiltonian. In all cases, our results devi- 
ated no more than 0.1% from the exact ground-state en- 
ergies, obtained by diagonalizing the Hamiltonian with 
the 0XBASH code [13]. For example for 29 Si we found a 
ground-state energy of —133.98 ± 0.04 MeV compared 
with the exact result of —133.95 MeV. Our method also 
reproduced correctly the ground-state spin in all cases. 

We applied our method to nuclei in the (pf+gg/2) shell, 
using the isospin-conserving Hamiltonian of Ref. 
Typical results are demonstrated in Fig. [TJ in which the 
absolute value of the logarithm of the Green's functions 
for the neutron orbital v — lpa/2 in 56 Fe (r > 0) and 
in 58 Fe (r < 0) are plotted versus |r| for /? = 4 MeV -1 . 
The linear fits (solid lines) were used in the calculation 
of the ground-state energy of 57 Fe. The deviations from 
the linear fits are shown in the insets of Fig. [TJ 

A direct application of the SMMC method to the odd- 
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FIG. 2: The energy of the 57 Fe nucleus calculated from the 
present method and direct SMMC are shown by solid and 
open squares, respectively. The error bars describe the statis- 
tical errors. Inset: the statistical errors for the energy of 57 Fe 
in the present method (solid squares) and in direct SMMC 
calculations (open squares) are shown on a logarithmic scale. 
The statistical errors for the energy of 56 Fe using the same 
Hamiltonian are shown by open circles. 



particle systems suffers from a sign problem which leads 
to very large statistical errors at low temperatures. In 
contrast, the method presented here does not have such 
problem. This is illustrated in Fig. [5] where we com- 
pare the energy and its statistical error for the 57 Fe nu- 
cleus in the present method (using the neutron Green's 
functions of 56 Fe) with the results obtained from the di- 
rect method. The errors in the present method remain 
roughly constant with /3. At j3 = 3 MeV -1 the statistical 
error in the direct method is about 5 times larger than 
the present method while at j3 = 4 MeV -1 it is about 20 
times larger. The inset shows the statistical errors on a 
logarithmic scale. For comparison we have also included 
the statistical error in the energy of the even-even nucleus 
56 Fe using the same Hamiltonian. 

We applied our Green's function method for fami- 
lies of odd-neutron isotopes: 47 " 49 Ti, 51-57 Cr, 53-61 Fe, 
59-65 Ni, 63-67 Zn and 71-73 Ge. The ground-state spins 
we determine are in agreement with experimental values 
in all cases except for 47 Ti, 57 Fe and 63 Ni. The anomalous 
ground-state spin of 57 Fe from the shell model perspec- 
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tive is well documented in the literature 

In our method we extract directly the odd-even ground 
state energy differences, and therefore this method is par- 
ticularly suitable for accurate calculations of pairing gaps 
(i.e., odd-even staggering of masses). 

When extracting an odd-even ground-state energy dif- 
ference such as AE m i n (A+) we use the Hamiltonian of the 
A+ nucleus for both the A+ and A nuclei. Since the fp+ 
gg /2-shell Hamiltonian we use is nucleus-dependent [5( , it 
is necessary to correct the ground-state energy of the A 
nucleus. As the latter is an even-even nucleus, this cor- 
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FIG. 3: Neutron pairing gaps A n as a function of mass num- 
ber A in fp + <?9/2-shell nuclei. The gaps calculated with the 
present Green's function method (solid circles connected by 
solid lines) are compared with the experimental gaps (open 
circles connected by dashed lines). The theoretical statistical 
errors are smaller than the size of the symbols. 



rection can be found in direct SMMC calculations for 
the A nucleus. However, this correction can also be esti- 
mated as follows. The dependence of the interaction on 
the nucleus is rather weak; the strengths of the multipolc- 
multipole interactions depend weakly on the mass num- 
ber A (oc A -1 / 3 ) and the monopole pairing strength is 
constant through the shell. The largest variation among 
neighboring nuclei is that of the single-particle energies 
£fj,(A) of the orbitals /i. Correcting for this variation, the 
neutron separation energy for the A+ nucleus is given by 



S n {A+) = -AE min (A+) + Y^[ e M ~ 



(5) 

where {n^A are the average occupation numbers for the 
A nucleus using the Hamiltonian for the A+ nucleus. The 
second term on the r.h.s. of ([5]) approximates the differ- 
ence between the ground-state energies of the A nucleus 
when calculated using the respective Hamiltonians for 
the A and A+ nuclei. We verified (in s<i-shell nuclei) that 
this approximation is highly accurate and well within a 
typical statistical error. In our calculations we used (O 
since the resulting statistical error is much smaller than 
the statistical error of direct SMMC calculations. 

The neutron separation energy for the A nucleus is 
given by a similar expression. The neutron pairing gaps 
can then be calculated from the differences of separation 
energies A n (A) = {-) N [S n (A+) - S n (A)}/2, where A 
can now be either an even-even or an odd-even nucleus. 

Our calculated pairing gaps are shown in Fig. [3] (solid 



circles) , where they are compared with the experimental 
pairing gaps (open circles) as determined from odd-even 
staggering of binding energies. Our results agree quite 
well with the experimental values; in most cases the de- 
viation of the theoretical pairing gaps from their experi- 
mental counterparts is less than 15%. Systematic devia- 
tions are observed for the iron isotopes above A — 59 and 
for the germanium isotopes. For the germanium isotopes 
the size of the model space might be insufficient, while 
the deviation for the iron isotopes indicates the necessity 
to refine our isospin-conserving Hamiltonian. 

Conclusion. We have described a practical method 
that circumvents a sign problem for calculating the 
ground-state energy of odd-particle systems in the shell 
model Monte Carlo approach. We have demonstrated 
the usefulness of the method by calculating pairing gaps 
of nuclei in the/p + 59/2 shell. This method can also 
be applied to other finite-size many body systems such 
as trapped cold atoms. In principle this method can be 
used more generally to calculate the lowest energy state 
for a given spin. However, when such a state is an ex- 
cited state, the statistical errors are larger and it is more 
difficult to identify the asymptotic regime. 
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